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Abstract 

The production of hadrons in heavy-ion collisions at RHIC in the low transverse-momentum 
(pt) region is investigated in the recombination model with emphasis on the effects of minijets on 
the azimuthal anisotropy. Since the study is mainly on the hadronization of partons at late time, 
the fluid picture is not used to trace the evolution of the system. The inclusive distributions at 
low pt are determined as the recombination products of thermal partons. The pt dependencies of 
both pion and proton have a common exponential factor apart from other dissimilar kinematic and 
resonance factors, because they are inherited from the same pool of thermal partons. Instead of 
the usual description based on hydrodynamics, the azimuthal anisotropy of the produced hadrons 
is explained as the consequence of the effects of minijets, either indirectly through the recombi- 
nation of enhanced thermal partons in the vicinity of the trajectories of the semihard partons, 
or directly through thermal-shower recombination. Although our investigation is focussed on the 
single-particle distribution at midrapidity, we give reasons why a component in that distribution 
can be identified with the ridge, which together with the second harmonic V2 is due to the semihard 
partons created near the medium surface that lead to calculable anisotropy in <f>. It is shown that 
the higher azimuthal harmonics, v n , can also be well reproduced without reference to flow. The pt 
and centrality dependencies of the higher harmonics are prescribed by the interplay between TT 
and TS recombination components. The implication of the success of this drastic departure from 
the conventional approach is discussed. 
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I. INTRODUCTION 



As the data on single-particle distributions of identified hadrons produced in heavy-ion 
collisions become more abundant and precise jl-7], more demands are put on theoretical 
models to reproduce them. It is generally recognized that in Au-Au collisions at y^s^T/v = 

ow transverse-momentum (p?) 



-Il2j. Because it has attained 



200 GeV at the Relativistic Heavy-Ion Collider (RHIC) the 
region (jpx < 2 GeV/c) is well described by hydrodynamics 
the status of the conventional approach, it is of interest to point out that an alternative 
approach without using the fluid description can also reproduce the same data on p T and 
azimuthal angular (0) dependencies and without using more parameters. The approach that 
we present here does not have the virtue of tracking the time evolution of the dense system, 
but it presents a different point of view on the origin of the azimuthal asymmetry. Instead of 
emphasizing early thermalization and the effects of anisotropic pressure gradient, we consider 
the non-flow effects of minijets due to semihard scattering of partons. The basic idea has been 



discussed previously [13l-ll5| . Here we present a more detailed phenomenological analysis of 



the pt and distributions of pions and protons produced at RHIC for a range of centralities. 

Our calculations are based on the recombination model [if]], which is a particular imple- 
mentation of the general approach of coalescence that has been shown to be successful in the 
intermediate region, 2 < p? < 6 GeV/c 17H19j| . Here we push to the lower region < 2 
GeV/c, in which the thermal partons are dominant. However, because semihard partons 
can lose energy to the thermal medium and result in local enhancement that is azimuthally 
anisotropic, there are non-trivial complications in the thermal sector. The point that mo- 
tivates our study is related to the question of what happens to the initial system within 1 
fm/c after collision. Semihard partons created within 1 fm from the surface will have already 
left the initial overlap region before thermalization is complete. They are the minijets that 
can give rise to dependence, not accounted for by conventional hydrodynamics. When the 
parton kx is low enough so that minijets are copiously produced, the corresponding effect on 
the anisotropy can become dominant, and is insensitive to the type of hadron produced. 

Another area of concern is the pt distributions in the low-pr region, where pion and 
proton appear empirically to have different behaviors. In the parton recombination model 
the hadrons should have the same inverse slope as that of the coalescing quarks if the hadrons 
are formed by recombination of the thermal partons, but because of the difference in the 
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meson and baryon wave functions, the net pt distributions turn out to be different. This line 
of analysis takes into account the quark degree of freedom just before hadronization, which 
is overlooked by the fluid description of the flow effect. The burden is to show that the data 
on t>2(pr) can be reproduced for both pion and proton at low p? in an approach using a 
common inverse slope T without relying explicitly on the hydro description of elliptic flow. 

The basic assertion in our study is that the recombination of thermal partons has two 
components, one is azimuthally isotropic, called the Base, while the other one is identified 
with the Ridge, which has the <fi dependence that is calculable. Our focus is on single-particle 
distributions at midrapidity and low p? and on how they are affected by minijets. The 
semihard partons that give rise to the observed minijets generate also the second component 
in the inclusive distribution. It will be our main task to show that the second component 
exhibits the properties of both the ridge and the second harmonic in <p. 

Recently there are experimental and theoretical studies of higher harmonic coefficients, 



v n , of the azimuthal asymmetry that 



rave been related to the fluctuations of the initial 



configuration of the collision system 2CH25l| . Although the phenomena are not of first order 



in importance compared to the second harmonics v 2, which has been regarded as the primary 
evidence for hydrodynamical flow [26], it is imperative for us to explain their origin in our 
approach that has no explicit formalism to connect the initial and final states. We shall 
show that their dependence on p? and centrality can be well accounted for by the thermal- 
shower component of recombination, so the minijets themselves cause the <fi anisotropy that 
leads to v n , whereas the dominant phenomenon in v<i is due to the enhancement in the 
thermal-thermal sector. 

We are aware that our approach is not in the mainstream and that we do not have a 
code to simulate the evolution of the dense system. However, it is of some value to have 
explicit analytic expressions that show why the pion and proton distributions in p? and 
(f) have common properties based on the parton distributions before hadronization, and 
how minijets can affect azimuthal harmonics in ways that are in accord with the measured 
behaviors. To have phenomenological evidence for the validity of an alternative approach 
that does not rely on hydrodynamical flow enriches the scope of inquiry into the various 
processes that can be important in heavy-ion collisions and may even cast doubt on the 
validity of the assumptions made in the conventional approach. 

Before entering into the details of our formalism, it is helpful to clarify possible confusions 
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of what we do with the conventional approach when common terms are used with different 
meanings. In hydrodynamical treatment of the dense system, it is usually assumed that rapid 
thermalization is completed in less than 1 fm/c and that the expansion of the system can 
then be described by the hydro equations with suitable assumptions about equations of state 
and viscosity Thus the words hydro and thermal are almost synonymous. In our treatment 
we use thermal without implying hydro. That is because we apply the notion of thermal 
to the soft partons only at late time just before hadronization. If the system takes over 5 
fm/c to equilibrate fully, that would invalidate the use of hydrodynamics from early time, 
but would not affect the validity of our approach. Or, if minijets make enhanced thermal 
contribution to the soft sector through energy loss of semihard partons to the medium 
without being a part of the equilibrated system from the beginning, then the thermal sector 
describable by hydro consists of only a portion of the soft hadrons in the final state, leaving 
another soft (but thermal) component that is outside hydro. These are possibilities that do 
not invalidate our approach; indeed, our treatment is aimed at accounting for the effects of 
those minijets. 



II. COMMON FORM OF HADRONIC SPECTRA 



We begin with a recapitulation of our description of single-particle distribution [16]. At 
low Pt we consider first the recombination of thermal partons, for which the pion and proton 
spectra at y = are given by 
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where T(qi) is the thermal distribution of the quark (or antiquark) with momentum qi, and 
IZh is the recombination function (RF) for h = it or p. On the assumption that collinear 
quarks make the dominant contribution to the coalescence process (so that the integrals are 
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one-dimensional for each quark along the direction of the hadron), the RFs are 



^7r(?i,?2,Pr) = I — - 1 ) , (3) 

Pt \~( Pt ' 



^p(?i,92,93,Pr) = / P (—,—,— ) 5 V — - 1 

VPt Pt Pt/ Pt 
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(4) 



where the details of fpiqi/pr) that depends on the proton wave function are given in 16], 
and need not be repeated here. The main point to be made here is that if the quark thermal 
distribution T(qi) has the canonical invariant form 

T(ft) = q^ = C qi e-^ T , (5) 
dqi 

where C has the dimension of inverse momentum, then the (^-functions in the RFs require 
that the hadron distributions p°dN h /dpx in Eqs. (1) and (2) have the common exponential 
factor, exp(— p^/T), for both h = n and p. The factors before the exponentials are different. 
The integrals in (1) and (2) must yield on dimensional grounds C 2 p\ and C 3 p T , respectively, 
apart from different multiplicative constants. Upon dividing the results of the integration 
by PoPt we get the general form 

dN TT(T) 

^V" =M h (pr)e-* r ' T , (6) 
p T dp T 

where, for y = 0, we set p = pr for pion and po = mr for proton, so that 

AC = N*C\ M p {p T ) = N?C 3 ^, m T = (p 2 T + m^l\ (7) 

TTLj' 

Nq and Nq being constants. Note that the factor p^/rriT in the proton spectrum causes 
the p/tt ratio to vanish as pt — > on the one hand, but to become large, as pt increases, 
on the other. When p? exceeds 2 GeV/c, shower partons become dominant and the above 
description must be corrected by the effects of thermal-shower recombination that limits the 
increase of the p/n ratio to a maximum of about 1 

Remaining in the low-pr region, p^ < 2 GeV/c, we want to demonstrate that a common 
value of T is shared by dN h / pxdpr for both h = it and p. The normalization factor M h {pT) 
in Eq. (jHJ) depends on centrality, which is a subject discussed in the Appendix. Here we 
consider a specific centrality, 20-30%, and fit the pt dependence of the proton spectrum 
using Eqs. ()6]) and ((7|) with free adjustment of the normalization constant. Figure 1 shows 
the result upon using 

T = 0.283 GeV. (8) 



The one-parameter fit (apart from normalization) is evidently very good compared to the 
data from Ref. It demonstrates that the proton is produced in that px range by thermal 
partons and that the flattening of the spectrum at low p T is due to the prefactor p^/rriT 
arising from proton recombination. 
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FIG. 1: Proton spectrum at y ~ averaged over <j> (hence, no 1/2tt factor) at 20-30% centrality. 
The solid line is a fit of the data by Eqs. © and ([7]) using T = 0.283 GeV with free adjustment of 
normalization. The data are from Ref. 

As already discussed at the end of the preceding section, the thermal parton distribution 
we consider is for the time just before hadronization. T in Eq. ([5]) is the inverse slope that 
we have determined here phenomenologically without the assumption that hydro description 
is appropriate for the entire period from collision to hadronization. It is, however, assumed 
that local equilibration is achieved for the soft sector at late time to justify the use of Eq. 
fl5]) for all pt < 2 GeV/c. We refer to T as inverse slope, instead of temperature, because we 
allow the possibility that the value of T can be affected by the motion of the collective system 
and become larger than the temperature defined in the local rest frame. For that reason 
the value of T should not be identified with what is referred to as freeze-out temperature in 
some fluid description. 

Having determined T, we have no more freedom to adjust the exponential behavior of 
the pion spectrum dN^/p^dpT- We show in Fig. 2 the data from PHENIX jl] on the pion 
distribution for 20-30% centrality; the exp(— pr/T) factor is shown by the solid line, the 
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normalization being adjusted to fit. The excellent agreement thus supports the assertion 
that both proton and pion spectra can be described by the same T in the exponential 
factor, exp(— px/T). For px < 1 GeV/c the pion spectrum is dominated by the resonance 
contribution which we cannot calculate for lack of knowledge about the RFs of hadrons 
above the ground states with orbital excitation. For that reason we show only the data in 
the region px > 1 GeV/c, which is sufficient to verify the commonality of T in the calculable 
part of our approach. 
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FIG. 2: The data for the pion spectrum are from Ref. The solid line that shows e PT / T has 
the same T as for proton. 



Our basic formulas for recombination shown in Eqs. ([TJ and (J2J) are valid for the direct 
production of all hadrons from thermal partons at any px- However, those hadrons are not 
the only ones produced in certain p T regions due to other processes, such as TS and SS 
recombination, which dominate at px > 2 GeV/c. For proton production TTT recombina- 
tion is prevalent at px < 2 GeV/c. For pions from resonance decays that are important for 
px < 1 GeV/c, Eq. (CQ) is inadequate to describe them. When the contributions from shower 
partons are significant even at px < 2 GeV/c, appropriate equations will be given below to 
describe them, px = 2 GeV/c is the upper limit of what we consider in this paper. 
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III. AZIMUTHAL ANISOTROPY WITHOUT HYDRODYNAMICS 



We now broaden our consideration to include dependence. For non-central collisions 
the almond-shaped average initial configuration leads to anisotropy. The conventional 
description in terms of hydrodynamics relates the momentum anisotropy to the variation of 



pressure gradient at early times upon equilibration [26|. The success in obtaining the large 
v 2 as observed gives credibility to that approach. Our alternative approach can be justified 
on the same basis that we can also reproduce the empirical V2, as we shall show. Concerns 
about higher harmonics are at a higher level of details which we shall also address. Our main 
objective is first and foremost to demonstrate that the essential phenomeno logical features of 
hadron production can be reproduced by considering the effects of semihard partons without 
explicitly treating the fluid flow. 

To include the azimuthal dependence, let us use p h (pT, 0, b) to denote the single-particle 
distribution of hadron h produced at mid-rapidity in heavy-ion collision at impact parameter 
b, i.e., 

where N part is the number of participants related to b in the known way through Glauber 
description of nuclear collision [27]. The main point of our approach is to assert that p h at 
low px can be separated into three components 

p h (p T , 0, b) = B h ( PT , b) + R h ( PT , 0, b) + M h (p T , 0, 6), (10) 

referred to as Base, Ridge and Minijet components, respectively. B h (pT,b) is azimuthally 
isotropic, while R h (pt, 0, b) and M h (pt, 0, b) are dependent. The first two components are 
due to the recombination of thermal partons (TT for pion and TTT for proton), while the 
third is due to thermal-shower recombination (TS and TTS) [16j. The latter is dominant in 
the intermediate pt region (2 < p T < 6 GeV/c), but is not negligible at low pt [151 ] . In a 
figure later on in this paper (Fig. 6) the relative magnitudes of the three terms will be shown. 
Because of the smallness of M h (pT, 0, b) relative to the other two for px < 2 GeV/c we shall 
ignore the shower partons in this and the following sections in order to emphasize the effect 
of semihard partons on the thermal sector. At the sacrifice of accuracy in reproducing t>2, 
the role of thermal partons in the ridge formation becomes clearer. More accurate result 
that includes the shower partons will be presented in Sec. VI. 
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Semihard partons created near the surface, and directed outward, can give rise to 
anisotropy in the thermal component. That is because each such semihard parton loses 
some energy to the medium, thereby enhancing the thermal motion of the soft partons 
near its trajectory Those thermal partons eventually lead to hadrons at late time that are 



dependent on the azimuthal angle of the semihard parton. In Refs. [13Hl5| those hadrons 
are identified with the ridge that stands above the background with characteristic peaking 
in 0, which is described by R h (pT,4>,b) in Eq. ffTOl) with improved treatment to be detailed 
below. The recoil semihard partons being directed inward are absorbed and randomized. 
They become a part of the medium consisting of all the soft and semihard partons that 
are farther away from the surface and are unable to lead to hadrons with distinctive 
dependence. That is described by the B h (px,b) component, which should not be confused 
with the 0-dependent bulk distribution in hydro description. For that reason we avoid using 
the term bulk. 

Ridge is a phenomenon characterized by an extended range in pseudorapidity rj and a 
narrow range in 28]. It may seem hard to relate the ridge to the term R h (pT, 0, b) in Eq. 
( ITUl) that has no r\ dependence. It is therefore important to emphasize that we are addressing 
here the dependence of the ridge at small rj, leaving aside the issues related to the range in 
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More discussions specifically 



7] that have been considered in our approach in Refs. 
on the ridge characteristics are given in Sec. V below. 

Earlier it was found that the azimuthal correlation between a semihard parton and a 
ridge hadron formed by TT recombination can be described by a Gaussian distribution in 
with a width a = 0.33 in order to reproduce the ridge data 3l|. That result was then 
extended in a study of the dependence of the ridge yield per trigger on the trigger angle 



relative to the reaction plane 



15 



321 ] . The key piece of physics that succinctly captures 



the essence of the correlation involving either a trigger or an untriggered semihard parton 
that generates the ridge is a quantity called 5(0,6). It is the segment on the initial ellipse 
through which semihard partons should be emitted if it is to contribute to the formation 
of any ridge particle that is directed at 0. The importance of 5(0, b) is that it relates the 
spatial and momentum asymmetries. The derivation of 5(0, b) given in 15] is based on 
the simple geometry of the average initial configuration taken to be an ellipse with width 
w and height h, where w — 1 — b/2 and h — (1 — b 2 /4) 1 / 2 in units of nuclear radius Ra- 
The fluctuations from that configuration will be considered later, but for now it is more 
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important to focus on the relationship between the spatial asymmetry (the short side of the 
ellipse being on the x axis) and momentum asymmetry of the emitted hadron ((p x ) > (Py))- 
Without repeating the derivation here, let us first state the result, and then follow up with 
a discussion on the physics involved. The result is 

S((f>,b) = h[E(e 2 ,a)-E(ex,a)], (11) 

where E(6i, a) is the elliptic integral of the second kind with a = 1 — w 2 /h 2 and 

6i = tan -1 ( — tan0j J , (pi = (p — cr, 2 = + 0", (12) 



for <pi < 7r/2, and an analytic continuation of it for (p 2 > tt/2. Thus S((p,b) is a compact 
formula that is completely calculable for any given b, and has the property that it decreases 
as <p is increased from to tt/2, more so at medium or large b than at small b. Since it will 
prescribe the (p dependence of R h (pT,(p,b), the average hadron momentum for non-central 
collisions is larger along p x than along p y . 

It should be noted that Eqs. (ITTj) and ffl2|) involve azimuthal angles only, which are 
meaningful in both the coordinate space and momentum space. The relationship between 
the angles is based on the correlation between the hadronic momentum in the ridge at angle 
(p and the direction <p s of the semihard parton emitted from the ellipse in the coordinate space 
with \(p — (p s \ < cr. Although the angles (p s and (p are at early and late times, respectively, 



they are nevertheless correlated, since a was determined in Ref. 



31j to fit the data on ridge 



yield as a function of the trigger angle <p s 32] . The only assumption here is that the trigger 



angle is identified with the semihard parton angle that is not directly measurable. Such an 
identification does not rely on the details of hydrodynamics. 

The physical origin of the <p anisotropy generated by semihard partons is that on the 
broader side of the spatial ellipse there can be more semihard partons within an angle a 
contributing to a hadron emitted with small <p at small y, where the curvature on the ellipse 
is small thus allowing a longer segment on the ellipse with normal in the range (p ± o. On 
the narrow side at the tip of the ellipse the curvature is larger thus restricting the segment 
through which semihard partons can contribute to a hadron at <p ~ tt/2. The mechanism 
that gives rise to this orthogonality between the spatial and momentum asymmetry axes is 
entirely different from that in the fluid description which is basically that the higher pressure 
gradient along the x axis in the initial state generates more hadronic momentum along that 
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direction in the final state. Without hydrodynamics we, of course, cannot describe the 
evolutionary history of the system. While the hydro approach assumes rapid thermalization, 
we allow unspecified time interval for expansion and equilibration except that by the time of 
hadronization at late time the soft partons have exponential px behavior, which is the only 
property we ascribe to the thermal partons, apart from the dependence of the enhanced 
thermal partons caused by the semihard partons on their way out of the medium at early 
time. If this approach can lead to sensible phenomenology of the azimuthal harmonics, as 
we shall show below, then it is an alternative that should be weighed against the merits of 
the conventional approach. 

The discussion above is about the <fi anisotropy of the thermal partons. Also to be 
considered is the role of the shower partons which are the fragmentation products of the 
semihard parton outside the medium before they hadronize. Inasmuch as the former reveals 
the effect of the semihard partons on the medium through which they traverse, the latter 
is the minijet manifestation of the semihard partons themselves by TS recombination. The 
hadrons formed are to be described by the third term M h (pT, <fi,b) in Eq. (fTUl) . The SS 



component is not considered in this paper because it is negligible at pr < 2 GeV/c 16 



331 ] . Shower partons can arise from semihard and hard partons created throughout the 

medium in random directions. Because of jet quenching the partons that emerge from 

the medium have reduced momenta, and the distribution of the shower partons generated 

by subsequent fragmentation peaks at low pt- They recombine with the thermal partons 

in the immediate vicinity of the emerging partons and therefore form hadrons that have 

approximately the same <fi angles as the initiating semihard or hard partons. Upon averaging 

over all events the azimuthal dependence of the TS term can have all harmonic components 

as in Fourier decomposition. Since the pt dependence of the azimuthal harmonics is what we 

shall calculate and compare with data, we summarize here the formulas for pion production 

by TS recombination that are relevant 

dNj s 2 f dpi dp 2 
p T dp T Pt 

where 



f ——T(pi)S(p 2 ^)n 7T (p 1 ,p 2 ,p T ), (13) 

J Pi P2 



S(P2,0 = J ^Y^Hq^SMq). (14) 

<S(p2) is the shower parton distribution integrated over all semihard parton q at the medium 
surface after momentum degradation parametrized by £, and Si{p 2 /q) is the distribution of 
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shower partons with momentum t>i i n a jet of type % with momentum q. The details of these 
quantities can be found in Refs. [la, |33 |. 



IV. SECOND HARMONIC OF ANISOTROPY 

This topic is usually referred to as elliptic flow, a terminology that is rooted in hydrody- 
namics. Since hydro is not the basis of our investigation, we use the more general language 
of harmonic analysis and refer to the harmonic coefficients 

.40*, 6) = <c 0S = C % c :Tf [P li b) ' < 15 > 

Jo «0P (Pt,0,o) 

where p h (pT,<f),b) in our formalism has the three components given in Eq. (IIP]) . We now 
describe the dependence of R h (pt, 0, £>) and M h (p T , 0, 6) separately. 

As discussed in the last section, R h (pj<,<f),b) contains the anisotropy arising from the 
initial elliptical spatial configuration through the 5(0, b) function that transforms the spatial 
to momentum asymmetry. We now insert some details omitted in our general discussion. 
Since the elliptical axes need not coincide with the reaction plane that contains the impact 
parameter vector b, we introduce a tilt angle ip2 and average over it. Furthermore, we modify 
the notation slightly by using 5*2(0,6) to denote what is defined in Eq. (fTTj) and write the 
average over 7r/2n as 

2 r /4 

-n/A 

We then define 5(0, fe) as the normalized 5 2 (0, 6), i.e. 



5 2 (0,&) = - f ' dfoSM-^b). (16) 

71" J-tt/4 



1 



2tt 



5(0, b) = £ 2 (0, b) / — / rf05 2 (0,6). (17) 



o 



Following our discussion in the last section on the ridge component of p h that responds to 
the minijets through TT recombination, we now can write 

R h ( PT ,<j ) ,b) = S( ( j ) ,b)R h (p T ,b), (18) 

where R h (pT,b) is the second of two components of dN^ T ^ /pxdpT- The exponential be- 
havior of the first component, which is the 0-independent base component B h (pT,b), has 
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a lower To than the overall T for the sum of the two thermal terms described by Eq. ([6]). 
Thus, with the unenhanced base thermal component expressed as 

B h (p T ,b)=Af h (pT,b)e^ To , (19) 

the enhanced ridge component is 

R h (p T , b) = Af h {p T , b)[e~^/ T - e-^ T °}. (20) 

We emphasize that the only factor that depends on the hadron type is Afh(pr,b). It is a 
specific property of the recombination model that the exponential factors of the hadrons 
(whether it or p) are inherited from those of the partons as discussed in the preceding 
section. Note also that T is the only unknown parameter introduced here. If for the present 
we neglect the TS component for the sake of simplicity, since it is small at low pt, we would 
have only the first two terms of p h (pT, 4>, b) in Eq. ([10]) . and the formalism up to this point 
should be sufficient to provide an approximate description of the second harmonic. 
Applying Eq. ( fl8l) to ( ITS"]) , we obtain for n = 2 

R h (p T ,b)± f*d<f> cos 205(0,6) 



v5(pT,b) 



B h (p T ,b) + R h (p T ,b) 
(cos 20) s 



where 



Z-\p T ) + l 



(cos20) 5 = J- / d<P cos 205(0, b), (22) 

Z{pT ) = ^ll = e PT/T> _ h T= JVI_ ( } 

yFT! B h (p T ) ' T-T y ! 

These equations are remarkable in that the b dependence resides entirely in Eq. (122]) and 
the pt dependence entirely in Eq. f|23|) ; furthermore, there is no explicit dependence on the 
hadron type. 

From Eqs. ( ITT]) . f|T6]) . f fTT]) and (|22]) we can calculate (cos20) 5 and obtain its dependence 
on b. The result is shown in Fig. 3(a). To check how realistic phenomenologically the fac- 
torizability of pt and b dependencies of v^ipr^b) is, we show in Fig. 3(b) the data from 
Ref. on t>2 (px, N paIt ) for three pt values, but shifted vertically so that they agree with 
the data for p T = 0.975 GeV/c for most of large iVp art . The diamond and square points are 



(21) 
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FIG. 3: (Color online) (a) Average of cos 20 weighted by S(cp, b) vs impact parameter b in units of 
Ra- (b) Common dependence of v^px, b) on N part for various px, shifted vertically for comparison. 
The diamond and square points are horizontally shifted slightly from the points in circles to aid 
visualization. The solid line is from (cos 2<f>) s shown in (a), but rescaled and plotted in terms of 
-^part- The data are from Ref. 2j. 

slightly shifted horizontally to spread out the overlapping points for the sake of visual distin- 
guishability. The fact that their dependencies on N part are so nearly identical is remarkable 
in itself. The solid line is a reproduction of the curve in Fig. 3(a) but plotted in terms of 
iVp art , and reduced in normalization by a factor 0.9 to facilitate the comparison with the 
data points. For iVp ar t > 100 the line agrees with the data on t>2 very well, thus proving 
the factorizability of Eq. (|2"T|) . For iVp art < 100, corresponding to &/-R4 > 1.3 or centrality 
> 40%, there is disagreement which is expected because the density is too low in peripheral 
collisions to justify the simple formula in Eq. ffTTj) that is based on no punch-through of 
recoil partons. 

Since T describes the pt dependence of the B h (p T , b) component, it is not directly related 
to any observable spectrum. Thus we turn to v%(pT,b) in Eq. ( l2Ti) for pion first and find 
that the low-pr data of v^ipr, b) can be well reproduced. In Fig. 4(a) is shown the data for 
pion from Ref. Q for 0-5% centrality; the solid line is the result of our calculation from Eq. 
d2U) using 

T = 0.245 GeV. (24) 
The fit though not perfect is remarkable because the normalization of is fixed by Eq. fl2T|) 
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without freedom of adjustment. Note that we have not used any more parameters besides T 
to accomplish this, which is a fitting procedure not more elaborate than the hydro approach 
where the initial condition and viscosity are adjusted. A better result will be shown below 
when the TS component is taken into account. 
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FIG. 4: v% at 0-5% centrality for (a) pion and (b) proton. The data are from Ref. 
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For proton v%(pt, b) we take the mass effect into account in the same way as we have 
done before [14j by working with transverse kinetic energy Ex, where 



E T (p T ) = m T (p T ) -m h , 
and adopt the ansatz that pt is to be replaced by E T in Eq. f[2Uj) . i.e., 

Z(p T ) = e^frrVT' _ L 



(25) 



(26) 



Using the same formula Eq. 02 ip without any change of T , we show in Fig. 4(b) that the 
result agrees with the proton data |2J fairly well. 

So far we have concentrated on v^ipr^b) at low px for 0-5% centrality as a first step 
toward demonstrating the phenomenological relevance of our approach based on thermal 
partons only in the first two terms of p h {pT, <P, b) in Eq. (|T0|) . To widen the Pt and b ranges, 
as well as to consider higher harmonics, we must include the third component generated by 
TS recombination. But before doing so, we pause in our study of the azimuthal anisotropy 
and revisit the subject of ridge. 
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V. THE RIDGE 



As remarked earlier in Sec. Ill, our study of the pt and dependencies of p h (pT, 0, b) at 
midrapidity in this paper does not give us the scope that can include the rapidity dependence. 
The Arj range of the ridge either in triggered correlation or in untriggered autocorrelation 



39- 



reaches as high as 4 34-381 and has therefore been regarded as lone-range correlation 

n n 

41] . That problem in the framework of our approach is addressed in Ref. [30J. The subject 
of our concern here is the property of the ridge at rj ~ 0; more specifically, we describe the 
effect of ridge in the inclusive distribution at low p?. 

Ignoring the third term in Eq. (I10p for the present discussion, and using Eq. (j!8p for the 
second term, we have 

p\p T) 0, b) = B h (p T , b) + 5(0, b)R h (p T , b) (27) 

so that upon averaging over 0, we obtain the two terms 

p h (p T , b) = B\p T , b) + R h (p T , b). (28) 

Their sum is the inclusive distribution with the exponential behavior given by Eq. (jBj) for 
Pt < 2 GeV/c, whereas they separately behave according to Eqs. (ITU1) and (12"U|) . B h (p T ,b) 
has been referred to as base, while R h (pt, b) describes the ridge. They are both the hadronic 
products of the recombination of thermal partons. 

It is not obvious by examining Eq. ( J20l) that R h (pT,b) exhibits ridge structure, but the 
derivation of 5(0, b) outline in Sec. Ill clearly indicates that R h (px, 0, b) has the quadrupole 
behavior because of the effect of semihard partons. That is, in addition to the 0-independent 
base term, the additional ridge term is made manifest at when semihard partons are within 
a cone of width a around 0, owing to the enhancement of the thermal partons in the cone due 
to energy loss by the semihard partons. The hadrons formed in the ridge has a higher (p^) 
than those in the base. For a single-particle distribution we have, of course, no trigger to 
select a direction around which the enhancement can be measured. But that does not mean 
that the effect of semihard partons (and therefore the ridge) is not present in the inclusive 
distribution. For b/RA = 1 we show in Fig. 5 the pr dependencies of B n (pT,b)/J\f n (pT,b) 
and R n (pT,b)/J\f n (pT,b), defined in Eqs. (fT9l and (1201) . by the (red) dash-dotted line and 
(blue) dashed line, respectively, and referred to as base and ridge. The former has T = 0.245 
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GeV according to Eq. ( T24"l) ; the latter is not exactly straight in Fig. 5 but can be fitted by 
exp(-p T /T R ) with 



T R = 0.32 GeV. 



(29) 



This larger inverse slope clearly indicates that the hadrons in the ridge are the products 
of enhanced thermal partons compared to those in the base. The sum p K {pT-,b)/M- K {pTib) ) 
which is the inclusive, is shown by the (black) solid line in Fig. 5, whose inverse slope is 
given by T in Eq. (|SJ). Note that Tr — T = 47 MeV is very close to the value 45 MeV 



that Putschke reportec 
inclusive distributions 



as the difference in the values of T between the triggered ridge and 
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FIG. 5: (Color online) The pt distributions of the base B 7T (pT,b)/M n (pT,b) in red dash-dotted 
line, ridge R n (pT, b)/M n (pT, b) in blue dashed line and their sum in black solid line. The centrality 
is for b = 1 in units of Ra- 

Figure 5 does not show the r\ and characteristics of the ridge, but the dependence of 
R n (pT, (ft, b) is completely contained in S(<f), b) as expressed in Eq. f|T8|) . Hence, ^(pt, b) and 
the ridge are intimately related, both being the consequences of semihard partons. If the 
shower partons generated by the semihard partons lead to a trigger particle, then the hadrons 
associated with that trigger would exhibit a peak in around the trigger direction, as was 



shown in Ref. 



311 ]. in agreement with data 32] . In fact, a prediction on the asymmetry of 



hadrons produced on the two sides of the trigger direction was subsequently verified to exist 
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in the data |42j. If the trigger direction is integrated over, then the distribution of the 
ridge hadrons would behave as R n (pT, <fi,b). 

VI. HIGHER HARMONICS 

In recent years the study of </> anisotropy has advanced from V2{pt, b) to higher harmonics 



43N46| . It is widely accepted that the coefficients v n with n > 2 are due to fluctuations of 



the initial configuration whose spatial eccentricity e n leads to the corresponding v n of the 
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It has therefore 



hadronic momentum distribution through hydrodynamical flow 
become our burden of proof that our approach can also reproduce the higher v n without 
flow. 

In the preceding section we have shown that v 2 can be understood in terms of the <fi 
dependence of the TT recombination of the thermal partons affected by the passage of 
semihard partons through the medium. Although we have no transport model to follow the 
evolution of the system, the space-momentum transformation is accomplished by studying 
the minijets emitted from the initial elliptical configuration, hence t>2. It is then natural 
for us to focus on the effects of the same minijets on the higher harmonics. In a sense 
minijets play a role similar to the fluctuations of the initial configuration, because their 
effects on the distribution present a departure from the consequence of the simple and 
smooth approximation of that configuration by ellipse, except that minijets are themselves 
fluctuations in the momentum space and do not depend on flow dynamics. In our formalism 
the minijets affect the low-pr region through TS recombination. 

Since minijets are produced in any given event in unpredictable directions, the average <fi 
distribution can have all terms in a harmonic analysis. The only aspect of the behavior that 
our formalism has a predictable power is the dependence on and centrality because the 4>- 
integrated TS component of recombination has already been formulated and parametrized. 
To be explicit, let us write the third component of p h {pr, </>, b) in Eq. ([TO]) as 

M h (p T ,0,b) = J(<P,b)M h (p T ,b), (30) 

where J(<j), b) describes the ^-dependent part of the minijet contribution, which is assumed 
to be factorizable from the average M h (pT, b) in the same manner as for R h (pT, 4>, b) in Eq. 
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([IS]) . Moreover, as in Eq. ( ITT)) . J((j),b) is the normalized form of J(4>,b) 

J(0, 6) = J(0, 6) y^i- jf #J(0, b) , (31) 

where J(0, b) contains all the harmonic components, cosn0, averaged over the tilt angle ip n , 
i.e., 

00 n r /2n 

J(0, 6) = 1 + b 2^ a n — / dip n cosn(0 — ^ n ). (32) 

„=2 ^ 



The b dependence in the above will be discussed below. We include n = 2 term in Eq. ( 132)) 
as an additional contribution to t>2 besides the one from Eq. ( ITS)) , which is dominant at low 
Pt because it is from TT recombination. With the M h term arising from TS recombination 
included, we shall go above the pt < 1 GeV region shown in Fig. 4. There is no way to 
calculate the amplitudes a n ,n > 2, but the px and b dependence of M h (p T ,b) is a unique 
attribute of our model, and will be put to test in our phenomenology of v n (pr, b) below. 
Including all three components of p h (pT, 4>, b) in Eq. ( fTUl) . we obtain from (|T5l) 



h f l\ _ (cos ncj)) s R h (px , b) + (cos n<fi) j M h (px , b) 
Vn[PT ' h) ~ p^b) ' (33) 

where 

p h (p T , b) = B h (p T , b) + R h {p T , b) + M h (p T , b), (34) 
1 f 27T 

(cos ncj)) j = — / d(fi cos ncj) J (4>, b). (35) 
2vr Jo 

(cosn<p) s is as defined in Eq. (|22|) for any n, but it is zero for n > 3 because of the periodicity 
of S(4>, b). Indeed, (cos nip) j receives contribution only from the a n term in Eq. ( 132]) because 
of the orthogonality of the harmonics. It is clear from Eq. (1331) that the pt and b dependencies 
of v^(pt, b) are no longer separable as in Eq. ( 12T|) . when M h {j>t, b) is included as is necessary 
for p T > 1 GeV. 

In the following we shall consider only pion production by TS recombination, since TTS 
and TSS recombination for proton is less important for p? < 2 GeV/c. The equation for 
dNj s /prdpr given in ( IT3l) and ( TL4l) are more elaborate than we need for M n (pT,b). The b 
dependence of £ given in [33[ is of a scaling form at intermediate pr and is complicated. We 
shall use the approximate form used in the earlier treatment [1G] where for the most central 
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collisions 



dNj s 



p T dpi 



rPr 

= — / d PlPl e-^ T S(p T - Pl ), (36) 
) Pt Jo 

POO 

S{P2) = ^effOg / dkkf g (k)S g (p 2 /k), (37) 

where only gluon jets are considered explicitly with f g (k) being the distribution of (semi)hard 
gluon created with momentum k. The factor a g = 1.2 is used to take into account the other 
(semi)hard partons whose contribution we approximate by adding 20% to the contribution 



from gluon jets [49J. The parameter £ e g is the effective fraction of (semi)hard partons 
created anywhere in the medium that emerges to fragment into clusters of shower partons; 
it is determined phenomenologically to be 0.07 [16]. For non-central collisions we regard 
M n (p T , b) to be proportional to C (N pSiTt ) N coi i(b) , where C(iV part ) oc N%£% [see Eqs. (Al) and 
(pjj)]; it is the normalization of the thermal parton distribution. N co \\(b) is the number of 
binary collisions that normalizes f g (k). We thus have 

C(iV part )iV coll (6) dN™ 



M*(p T ,b) 



(38) 

6=0 



C(N™?)N coll (b = 0)p T dp T 
which completely specifies the px dependence. There is just one more piece of physics that 
needs to be added. That is the decrease of average path length in the medium as the collision 
becomes more peripheral. Its consequence is that more fraction of the (semi)hard partons 



can emerge from the medium as b increases. In [33] a detailed study of dependence of the 
nuclear modification factor on and b has been carried out. For the purpose of promoting 
a transparent connection between the harmonic coefficients v n and the input, we make the 
simple first-order approximation here that the increase of minijets with b can be expressed 
as a linear rise shown in Eq. ( 13"2"|) . which then exhibits very plainly the <fr and b dependencies 
of J((j),b), and therefore also M 71 " (j>t, 4>,b). 

It is useful to have a visual comparison of the various components of p n {pT, b) in Eq. (|34|) . 
In Fig. 6 we show B n (pT, b), R n (pT, b) and M n (pT, b) for 0-10% centrality, determined from 
using Eqs. (fl9]l . ( 1201) and ( )38l) . respectively. Evidently, R n {j>T,b) and M n (pT,b) become 
increasingly more important at increasing px- They set the scale of v^(pT,b) through their 
roles in Eq. (13"3"]) . For specific harmonics, we limit ourselves to n = 2, 3 and 4 and calculate 
v 2(PT,b),v^(pT,b) and v%(pt, b). The results are shown by the solid lines in Fig. 7, where 
the data are from PHENIX 50j. The values of the parameters used are 

a 2 = 0.6, a 3 = 1.6, a 4 = 1.2. (39) 
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FIG. 6: (Color online) B" T {p T ,b), R n (p T ,b) and M n (p T ,b) for 0-10% centrality. 

It is remarkable how well the calculated curves agree with the data in p? dependence for 
four centrality bins in each case. One parameter a n for each n can affect only the magnitude 
of v n (pT,b), so the excellent reproduction of the pt and b dependencies reveals the basic 
attributes of the approach that we have taken to describe the harmonics. 

We note that Fig. 7(a) represents an improvement over Fig. 5(a) owing to our use of Eq. 
(|33|) instead of (I2T!) for v% (pt, b). It is clear that the difference is due to the inclusion of the 
third term in Eq. ( TlUj) . However, it is also significant to point out that the change is not 
large except in the ranges of pt and b. That is, the major features of v%(pt, b) can be well 
reproduced by considering TT recombination only with the neglect of the TS component. 

The fit of the data is not perfect for v^(pT,b) and v\(pT-,b) at 30-40% centrality. That 
may be due to larger fluctuations of the minijet contribution at less central collisions, for 
which more detailed study is called for. 

Note that the curvatures of the lines and the data for V2 are different from those of v% 
and v^. To see the origin of that difference, we rewrite Eq. ( 1331) to reflect the dominance of 
the numerator of v%(pr, b) by ^(pr, b), and of v 3 4 (pr, b) by M n (pT, b) 

R n (PT,b) 



v 2(PT,b) ~ (cos 20) 
<(PT,b) 



{cosm 



5 ^(Pr,6)' 
M*{p T ,b) 



n = 3,4. 



(40) 
(41) 



3 P n (PT,b) ' 

Thus the pt dependencies of v^ipr, b) and v% ±{pt, b) are dictated by those of 
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FIG. 7: pt dependencies of (a) v^ipr, b), (b) v^(px, b) and (c) vJ(pt, b) for four centralities in each 
case. Data are from Ref. [5o|. Solid lines are the results of our calculation using a n given in Eq. 

HMD- 
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R n (pT,b)/ p n (pT,b) and M n (pT,b)/p n (pT,b), respectively, which are prescribed by the be- 
haviors shown in Fig. 6. is convex upward because, as pt increases, M(px, b) becomes 
larger so the increase of R n (j)T, b)/p n (pT, b) decelerates; in fact, it would decrease as px gets 
above px ~ 3 GeV/c, a property that has previously been obtained in Ref. 14| because 
of TS dominance. On the other hand, v£ and v\ are concave upward because R n (pT,b) is 
much larger than M n (pT,b) around pr ~ 1 GeV/c so M n (j)T,b)/ p n (pT,b) is suppressed at 
low pt- Eventually, as pt gets much larger, SS term will become important and turn v J 4 
over and diminish them. Since the properties of the three components in Fig. 6 are specific 
results of the recombination model, the appropriate curvatures of the solid lines in Fig. 7 in 
agreement with the data lend support to our minijet approach to the treatment of azimuthal 
anisotropy. 

Our study here is mainly a demonstration of principle in that minijets are important and 
can explain all the low-p^ data in the recombination framework. However, it is important 
to bear in mind that what we have shown is the sufficiency of our approach to reproduce 
the data, but not necessity. Neither is the hydro approach necessary. Indeed, there is no 
theoretical treatment that can prove necessity. Nevertheless, it is significant to recognize that 
various dynamical mechanisms can be responsible for the same phenomenological features 
of the hadronic observables. By the same token, a combination of those mechanisms may 
be at play in reality. The base component in our description could possibly be treated 
by hydrodynamics if thermalization is rapid for the subsystem that is left behind after 
the emission of semihard parton near surface. There are other related issues concerning 
the fluctuation of the initial configuration and the variation of the thermalization time for 
various different eccentricities. Such complications combined with the effects of minijets that 
we have found here open up a range of possibilities, on which our present treatment may 
reveal only a restricted view that is opposite to the traditional hydro view. A comprehensive 
study that includes both components of the mixture is a worthwhile problem for the future. 
For now, our simple remark is that the common usage of the term "elliptic flow" for t>2 is 
inadequate in generality and tends to be misleading. 
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VII. CONCLUSION 



We have shown that the major properties of pion and proton production at low px can be 
reproduced in our formulation of hadronization that includes the effects of minijets. The pt 
spectra have exponential behavior, exp(— px/T), with a common value of T for both tt and p 
that is the same as the T of the thermal partons just before hadronization. Minijets generate 
azimuthal anisotropy both through energy loss to the medium and in creating shower partons 
that recombine with the thermal partons. Harmonic analysis of the ^-dependence leads to 
v n (pT,b) that agrees with the data. 

We have also shown that the ridge phenomenon is a consequence of minijets. Although 
our study in this paper is limited to the small-// region of inclusive distribution, the ridge 
component is shown to have a harder pt spectrum because of the enhancement of the thermal 
partons. The dependence around a trigger was described in [3l|, but now we show that 
when integrated over the trigger direction the ridge component in the inclusive distribution 
generates the quadrupole v 2 (pT,b) with the correct pt dependence. Thus, v 2 and the ridge 
are tightly related. 

Since our treatment is only for the system at late time, we employ no model to carry out 
the development of the system from early time. The thermal partons are assumed to have 
an exponential form that is determined by phenomenology. Not following the evolution of 
the system is not equivalent to an assumption that the system does not expand. It is just 
that we do not make any assumption concerning the equilibration time or the properties 
of the fluid. Obviously, we do assume that by the time of hadronization there is local 
thermalization to justify the use of T. The claim we make is that taking the minijets into 
account is sufficient to reproduce the measured azimuthal anisotropy. We cannot exclude 
the validity of hydrodynamical flow, but we do show that the phenomenology that supports 
the flow dynamics provides the same support for our approach. Thus there are two possible 
descriptions of the low-pr process, neither of which can claim exclusive validity. The reality 
may even be a combination of both. 

While further investigation is needed to determine the extent of the admixture of flow 
and minijets at RHIC, it is conceivable that in collisions at LHC the density of semihard 
partons is so high initially that the system has insufficient time for equilibration before the 
abundant minijets created near the surface dominate the expansion characteristics, even 
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though the higher density of soft gluons speeds up the thermalization process at the core of 
the medium. If that is so, then one may think of what we have done here as the lower-energy 
precursor of what needs to be done at higher energies. The study of the px spectra at LHC 



49]. 



has already shown the importance of minijets through thermal- shower recombination 
It will therefore be natural to apply the formalism developed here to elucidate the problem 
of azimuthal harmonics measured at LHC. 
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Appendix A: Centrality Dependence of Hadronic px Distributions 

Having obtained the correct centrality dependence of " in Fig. 5 of Sec. IV, which is 
totally calculable without free parameters, we consider here the centrality dependence of 
the inclusive spectra p h {pT, b). We note that the unknown normalization factor M h in Eqs. 
06]) never enter into the calculation of " because of cancellation, but for p h (pr, b) they must 
be reckoned with. As indicated in Eq. (J7J), A/" 77 and M p are proportional to C 2 and C 3 , 
respectively, due to qq and qqq recombination. The magnitude C of the thermal partons 
depends on in a way that cannot be reliably calculated By phenomenology on the pion 
spectrum it was previously estimated for p? > 1.2 GeV/c [15|], but that is inadequate for our 
purpose here; moreover, Af n and M p have different statistical factors that can depend on b 
because of resonances. We give here direct parametrizations of the normalization factors in 
Eq. ([7]) in terms of A part 

N-C 2 = 0.667JVJS;, (Al) 
A£C 3 = O.UW^l (A2) 

Using them in Eqs. ([6]) and (JTj) we obtain the curves in Fig. 8 (a) pion and (b) proton for 
three centrality bins. They agree with the data from PHENIX [l(] very well over the range 
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of pt shown. In all those curves T is kept fixed at 0.283 GeV, thus reaffirming our point 
that both pions and protons are produced by the same set of thermal partons despite the 
apparent differences in the shapes of their p T dependencies. 
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